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ON THE MAGNITUDES OF COMPACT SETS IN EUCLIDEAN 

SPACES 

JUAN ANTONIO BARCELO AND ANTHONY CARBERY 


Abstract. The notion of the magnitude of a metric space was introduced 
by Leinster in \m and developed in ca, US], E], [13 and m, but the 
magnitudes of familiar sets in Euclidean space are only understood in relatively 
few cases. In this paper we study the magnitudes of compact sets in Euclidean 
spaces. We first describe the asymptotics of the magnitude of such sets in both 
the small and large-scale regimes. We then consider the magnitudes of compact 
convex sets with nonempty interior in Euclidean spaces of odd dimension, 
and relate them to the boundary behaviour of solutions to certain naturally 
associated higher order elliptic boundary value problems in exterior domains. 
We carry out calculations leading to an algorithm for explicit evaluation of the 
magnitudes of balls, and this establishes the convex magnitude conjecture of 
Leinster and Willerton m in the special case of balls in dimension three. In 
general we show that the magnitude of an odd-dimensional ball is a rational 
function of its radius. In addition to Fourier-analytic and PDE techniques, the 
arguments also involve some combinatorial considerations. 


1. Introduction 

Motivated by considerations of a category-theoretic nature, Leinster m has intro¬ 
duced the notion of the magnitude of a metric space. Magnitude is an important 
new numerical invariant of a metric space which shares some of the more abstract 
properties of the Euler characteristic of a category (or of a topological space), and 
indeed both can be seen as special cases of the notion of the Euler characteristic or 
magnitude of an enriched category. In particular, the inclusion-exclusion principle 
enjoyed by the Euler characteristic provides important motivation for the hoped- 
for properties of magnitude. More generally, magnitude is designed to capture the 
“essential size” of a metric space in a more subtle way than cruder measures such as 
cardinality or diameter, and at the same time it will also contain further significant 
geometric information concerning the space. For a much more detailed discussion 
of these issues see E], ffu, m and [T5] . 

Leinster’s definition of the magnitude of a finite metric space bears close resem¬ 
blance to notions of a potential-theoretic nature, and Meckes [16] and m has de¬ 
veloped this perspective to the point where a tractable definition of the magnitude 
of a positive-definite compact metric space can now be given in terms analogous to 
those of classical capacity. This provides the starting point for our investigations. 

Before describing our results, we give a little more informal background on magni¬ 
tude in order that our contributions can be placed in context. 
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1.1. Definitions of magnitude and connection with capacity. Given a finite 
metric space {X, d ), Leinster [11] defined its magnitude as the value 

1^1 = w{x) 

xeX 

whenever w : X —>■ R satisfies 


for all X € X. 


Y = 1 

yeA 


It is easy to check that any two such w will give the same value for the magnitude, 
and if no such w exists we declare the magnitude to be undefined. Under the 
mild additional condition that X be positive-definite (meaning that the matrix 
ig positive-definite), its magnitude is defined, see [11], [16] and [T7] . 
Finite subsets of Euclidean spaces are always positive-definite, m- It is easy to 
check (see below for the argument) that |0| = 0, |{t}| = 1 and that if Xjy is an 
iV-point space consisting of the vertices of a simplex in R" (with n > iV — 1) which 
are all equidistant t from one another, then 


( 1 ) 




N 

1-f (IV - l)e-‘■ 


We describe a compact metric space as positive-definite if every finite subset is 
positive-definite, and the magnitude of a compact positive-definite metric space 
{X, d) has been defined [16] as 

\X\ = sup{|S| : S a finite subset of X}. 

(Once again, any compact subset of a Euclidean space is positive-definite.) 

Bearing in mind Leinster’s definition, it is natural to consider the class of finite 
signed Borel measures /r on a compact positive-definite metric space X, and for 
such a p, introduce its potential function given by 

d>^(a:) := [ 

Jx 

If there is such a finite signed Borel measure fj, on X satisfying 

$^(x)= [ e-^^^'y'>dy{y) = lonX, 

Jx 

then y is called a weight measure for X. It is known ([I6],[I2]) that if X is compact 
and positive-definite, then |X| = y{,X) for any weight measure y. 

If we take X = [—R, i?] C R with the usual metric, then one simply checks using 
integration by parts that, with dec denoting Lebesgue measure on R, + 

der|[-RR]) is a weight measure for X, and hence |[—= i?-|- 1. (See [TT] and 
m) This is the only example of a nontrivial compact convex set in a Euclidean 
space whose magnitude was hitherto known. 

In nearly all examples for which the magnitude is known explicitly one can fairly 
easily identify a weight measure; in particular when the metric space X enjoys a 
lot of symmetry one expects a weight measure to reflect this symmetry and this 
leads to a limited range of possibilities. For example, in the case of the simplex 
mentioned above, we just set w{y) = a for all y and observe that the equation 



MAGNITUDES OF COMPACT SETS 


3 


= 1 becomes a + {N — l)e *a = 1, leading to ([T]). For more 
examples see m. m, |20) . Our approach here is different and is motivated by 
connections with differential equations. 

There is a clear analogy between magnitude and the notion of capacity as developed 
in classical potential theory, and which we now describe in a very informal manner. 
A possible definition of the a-capacity of a compact metric space is 

cap„(A) = sup/r(A) 

where the sup is taken over all positive finite Borel measures y on X such that 

[ d{x,y)~°dy{y) < 1 on A. 

Jx 

The study of the potentials d{x,y)~°‘dfi{y) and the associated a-capacities has 
a long and distinguished history. When we are in Euclidean space M" with the 
usual metric, we can take advantage of the Fourier transform to characterise the 
a-capacity of a compact subset A for 0 < a < n as 

caPa(A) = c„,„inf | / |Cr"“|/(C)Pd^ : / > 1 on A 

where Cn^a is a certain dimensional constant and ^ denotes the Fourier transform. 
See for example [I] for a thorough and detailed discussion of potential theory and 
capacity in the Euclidean setting. 

Meckes in m develops the analogy between magnitude and capacity beyond the 
formal level, and gives an extremal characterisation of magnitude in some generality. 
In the case of compact sets A in Euclidean space K." this characterisation can be 
realised as 

(2) |A| = ^ inf {: / G / = 1 on a} 

where ujn is the volume of the unit ball in K” and i?™(]R”) is the Sobolev space 
of functions whose derivatives of order up to m are in L^(]R"). More precisely, 
is the space of Bessel potentials of order m, and its norm is given by 

II/IIIf™(B") ll(d ~ ^)™^^/IIl2(R")- 

(There are many different - but equivalent ~ norms on the space i7"*(]R") which 
are regularly employed in the literature, especially when m is an integer. We 
emphasise that throughout this paper, we exclusively use the definition above.) 
Note that when to = (n + l)/2 > n/2, functions in are continuous by the 
Sobolev embedding theorem so the prescription / = 1 on A makes sense pointwise. 

One can informally motivate this result of Meckes in line with classical potential 
theory as follows. Adopting the standard convention from Fourier analysis - which 
differs from that used by Meckes in [16] and m - that the Fourier transform is 
given by 

m = [ /(x)e-2--?da;, 

then it defines an isometry on and we have 

(3) e-l-l(^) = nlwn (l+ 471^1^1^) 


. + 1 
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(See for example Stein, [18].) So if /r is a weight measure for X C R” we can extend 
it to be a measure on the whole of K" in the canonical way, <l>;_j becomes defined 
on all of M" and we have d)^(a;) = 1 on X. On the other hand, by taking Fourier 
transforms, 

= n!a;„ (l + ^ /i(0- 

Therefore 

/j(X) = [ = [ [ e-^^-y^dn{y)dfi{x) 

JX 

= [ n!a;„ (l + dTT^I^p) ^ 

JR^ 

= ^[ lC(OP(l + 47r2|?P)'^d^ 

n'.UJn jRn 

and this last expression is exactly ;n^||4>^||^(„+i)/2(Kn)- 

It is easy to see using Plancherel’s theorem and the binomial theorem that when m 
is an integer the inner product (•, corresponding to the H'^ norm is given 

by 

(4) {/,ff)//™(R") := ^ D^f ■ Dig 

where f = for j even and f = VAb“i)/2j j:pj. j shall make 

systematic use of this in what follows. 

One should note that although the right-hand side of equation ([2]) apparently de¬ 
pends upon the Euclidean space in which X sits, the original definition of the 
magnitude of X is an intrinsic metric invariant, and so formula ([2]) will hold for 
whatever Euclidean space M" in which we can embed X. 


1.2. Magnitude and geometric invariants. It is perhaps not a priori clear why, 
despite its ubiquity, the exponential function s i—>■ e“® is singled out to appear in the 
definition of magnitude in composition with the metric. Ultimately the reason for 
this is that its key property of converting addition to multiplication distinguishes 
it as the (essentially) unique function such that, if di{d{x,y)) is the kernel of the 
potential function, then the magnitude of a finite metric space (as dehned via ik) 
and the Euler characteristic of a finite category share a common generalisation in 
the Euler characteristic or magnitude of an enriched category. In particular, with 
'I'(s) = e“®, the triangle inequality gives us 

d!{d{x,y)) > 'd’{d{x,z))d!{d{z,y)) 

for all X, y, z e A; this inequality is to be thought of as analogous to the category- 
theoretic composition of a map from x to z with a map from z to y to obtain a 
map from x to y. For more details see m- 

In the category of finite sets, magnitude is just cardinality, which obviously satisfies 
the familiar inclusion-exclusion principle 


#{A \JB) + #{A n U) = #(A) + #{B). 
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Similarly, in the category of finite simplicial complexes, the classical Euler charac¬ 
teristic X can be characterised by x(0) = 0, x({2;}) = 1 and 

xiA UB)+ xiA r\B) = xiA) + x{B). 

However, magnitude in the setting of finite metric spaces is manifestly not finitely 
additive as the example of the simplex with equidistant vertices, o, indicates. Nev¬ 
ertheless, formula © tells us that as t —>■ oo, the magnitude of Xjv does asymptot¬ 
ically satisfy the inclusion-exclusion principle since |^jv| approaches N as t —>■ oo. 
More generally if (X, d) is any finite metric space and if we let tX = (X, td) for 
t > 0, it is known (see [TT]) that |tX| —>■ as t —>■ oo, and so inclusion-exclusion 
holds in an asymptotic sense for the class of finite metric spaces. On the other hand, 
if we consider the (admittedly somewhat limited) class /C of compact convex sets in 
K, then ii A, B G 1C are such that AU B G 1C, we have the exact inclusion-exclusion 
identity 

\AUB\ + \AnB\ = |H| + \B\ 

for magnitude (since the compact interval [a, b] has magnitude 1 + (b — a)/2 as we 
have seen above). These considerations lead one to hope that magnitude in general 
will satisfy the inclusion-exclusion principle in some asymptotic sense, and perhaps 
exactly so under certain more restricted circumstances - such as in the presence 
of convexity. (One may remark that the classical theory of capacity of compact 
convex sets is richer than the general theory, and from the category-theoretic point 
of view convex sets are natural in so far as the path of minimal distance between 
any two points is contained in the set.) 

Examples of metric spaces in which there is a lot of symmetry, and therefore where 
it is possible to evaluate the magnitude directly, yield further insight into some of 
the geometric characteristics which magnitude might capture. In analogy with the 
calculation for the simplex above, if X is a compact metric space and C? is a compact 
group with Haar measure u, acting transitively and isometrically on X, then the 
only sensible candidate for a weight measure for tX = (X, td) should have constant 
density \a{t) (for some A > 0) where a{t) = 1 (the integral is 

independent of x and y). So the magnitude of tX should be 

I I Jge-*'^G’92')du(g) 

Exact evaluation or asymptotic analysis of such expressions as t ^ oo has led to 
the following results: 

Theorem A (Willerton, [20). 1 Let X be a two-dimensional homogeneous compact 
manifold. Then 

\tX\ = ■^Aiea{X)t‘^ x(X) -b 0(t"^) 

Ztt 

as t ^ oo, where x is the classical Euler characteristic. 

Theorem B (Willerton, [20). 1 Let n > 1. Then there is an explicit formula for 
|tS”| with its geodesic metric, which yields 

= _^Vol„(S")r -b -^tsc(S’^)t”-2 -b • • • -b x(§”) + 0(e-‘) 
nluJn GnlUJn 

as t ^ oo, where Vol„ denotes n-dimensional volume, x is ths classical Euler 
characteristic and tsc is the total scalar curvature. 
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(Note that x(§”) is 2 for n even and 0 for n odd.) Willerton also has formulae for 
the magnitude of spheres in Euclidean spaces with the subspace metric which are 
explicit when n = 1 and 2 and give the first few terms in the asymptotic expansion 
as t —> oo for larger values of n. In addition, he has obtained the first four terms 
in the asymptotic expansion for a general n-dimensional homogeneous compact 
manifold. 

Together with the result that for a nonempty compact interval / C K we have 

|t/| = iLength(/)t + x(/) 

(where x{I)j the Euler characteristic of I, equals 1), these results of Willerton sug¬ 
gest that 1 1 —>■ \tX\ should have a polynomial flavour as t —)• oo, with the coefficients 
of the polynomial part encapsulating important geometric information concerning 
X. Combined with the results of numerical calculations [TO] and some parallel re¬ 
sults when the ^^-metric is used instead of the usual metric (see m), they provide 
evidence for the following conjectures of Leinster and Willerton from [12] : 

Conjecture 1 (Leinster-Willerton). For suitable (i.e. so that perimeter and clas¬ 
sical Euler characteristic are at least defined) compact X C we have, as t ^ oo, 

\tX\ = ;^Area(X)t^ -I- iperim(A)t -|- x{^) + o(l)- 


Here we see that intermediate-dimensional geometrical characteristics (in this case 
the perimeter) of X are expected to play a role. We shall not be directly concerned 
with Conjecture |T] in the current paper, but will instead focus on the following, 
known as the convex magnitude conjecture of Leinster and Willerton, m- 


Conjecture 2 (Leinster-Willerton). Suppose X C K" is compact and convex. 
Then 1 1 —> \tX\ is a polynomial of degree n and moreover 


(5) \tX\ = 

n\uJn 


Surf(aX) 

2(n - l)!a;n_i 


± 71—1 


■ • +1 = E — 

^ llLVi 


i=0 


where Vi{X) is the i ’th intrinsic volume of X and uJi is the volume of the unit ball 
in RL 


The intrinsic volumes V) are classical integral-geometric invariants defined on the 
class of compact convex sets K. in Euclidean space, and Vi{X) encapsulates the 
quantitative f-dimensional information concerning X. In particular we have Vi{tX) = 
fVi{X). Indeed, Vn is ordinary volume, Idi-i is half the surface area, and Vi for 
I < j < n — 2 captures the f-dimensional information on the boundary dX, while 
Vo (A) = 1 for nonempty X and Vo{X) = 0 when X is empty. Each V) is a valuation 
on 1C, which means that Vi(0) = 0, and ii A, B € 1C are such that AU B € 1C, then 
we have the inclusion-exclusion identity 

ViA UB) + V(A nB) = V(A) + V(B). 

For more detail see m- So the convex magnitude conjecture implies in particular 
that magnitude is a valuation on 1C. On the other hand, it is known by Hadwiger’s 
theorem (see again [TO]) that any valuation on K, which is continuous with respect to 
the Hausdorff metric and which is invariant under Euclidean motions must be given 
by a linear combination of the intrinsic volumes. It is also known (see m,m) that 
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a valuation which is translation-invariant and monotone is necessarily continuous. 
Clearly magnitude is invariant under Euclidean motions and is monotone, so the 
convex magnitude conjecture is true if and only if magnitude is a valuation and the 
conjecture is true for a suitable family of compact convex sets - such as the closed 
balls - which serves to normalise the precise coefficients of the polynomial arising. 
It is therefore of considerable interest to determine the validity of the conjecture in 
the special case of closed balls in R", and this issue is the main focus of the present 
paper. 


On the other hand, one does not see the inclusion-exclusion principle arising in 
classical potential theory, so perhaps a certain amount of caution is in order in 
approaching this issue. 


To fix ideas, the convex magnitude conjecture predicts that the magnitude of the 
closed ball of radius R in R" will be: 

1 


n = 1 
n = 2 


n = 3 


n = 4 


n = 5 


n = 6 


n = 7 


2 ! 


ttR 


r3 

— + R^ + 2R+1 


+ 


7ri?3 


i?5 i?4 

R^ ttR^ 

TT ^ 


3i?2 

2i?3 

5i?4 


3ttR 


4 

72 


-f 1 


57ri?3 

4ij3 


+ 1 
5i?2 


3i?^ 


15ttR 

16i? 


-HI 


etc. 


1.2.1. What is currently known about the convex magnitude conjecture. It is clear 
from © that 1 1 —> \tX\ is an increasing function of t when X is compact and convex. 
Meckes [T7] has shown that this function is continuous on (0, cxd), and Leinster [TT] 
(see also Meckes [H]) has shown that for arbitrary compact sets in R" we have 


1^1 > 


Vol(X) 

nluJn 


When n = 1 the conjecture is true (see above), but when n > 2, methods based 
entirely on symmetry - and indeed on weight measures - will not suffice to resolve 
it, and instead we turn to techniques of differential equations. 

We remark that a natural analogue of the convex magnitude conjecture does 
hold for convex bodies in R” when instead of the elP norm we use the norm. 
See Theorem 4.6(2) of [13]. 


1.3. Main results of the present paper. We shall first establish that the asymp¬ 
totics as t —>■ 0 and t ^ oo which are predicted by the Leinster-Willerton convex 
magnitude conjecture do indeed hold, and in fact do so more generally for nonempty 
compact sets. That is, we prove: 
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Theorem 1. Let X be a nonempty compact set in R”. Then 

I RX I —)• 1 as R ^ 0 

and 

R-^\RX\ asR-^oo. 

n'.UJn 

As a consequence of the first statement we have that t i-A \tX\ is continuous also at 
t = 0. 

We shall next be concerned with the validity of the convex magnitude conjecture 
for the class of convex bodies in Euclidean space K", that is, compact convex sets 
with nonempty interior. Our starting point will be Meckes’ formula ([2]) for the 
magnitude of such a convex body and we shall try to relate this to the formula ([S]) 
appearing in the convex magnitude conjecture. 

Meckes has already observed m that the extremiser in ([2]) exists and is unique, 
and also satisfies the associated Euler-Lagrange equation 

(6) (/ - A)("+i)/ 2/ = 0 weakly on R" \ X 

which, when n is odd, is an elliptic differential (as opposed to pseudodifferential) 
equation. In order to be able to work with differential rather than pseudodifferential 
equations we shall from Section^ onwards assume that n is odd. We shall first show 
that there is a unique member of which is a solution to ([S]) and which 

is identically 1 on X. This solution therefore extremises the energy ||/||^(„+i)/ 2 (][{n) 

over all / e with / = 1 on X. 

The convex magnitude conjecture predicts a single term corresponding to the vol¬ 
ume of X, with the remaining terms mainly relating to the boundary dX. Our 
next task is therefore to develop a formula for the extremal energy expressed in 
terms of Vol(X) and the values of the extremising function near dK. This is car¬ 
ried out in order to highlight the contribution to the extremal energy arising from 
the boundary of X and hence ultimately to facilitate a comparison between m and 
®. The formula is analogous to classical representations of capacities as boundary 
integrals of functions of solutions to the associated partial differential equations, 
and it appears in Theorem [5] below. 

We then turn to the problem of explicit identification of the extremiser in the case 
that X is a nontrivial closed ball. This is achieved by using spherical symmetry to 
reduce matters to consideration of certain ordinary differential equations. 

Finally, having obtained the extremiser, we use our formula from Theorem [5] to 
obtain an explicit evaluation of the magnitude of a closed ball. (Strictly speaking 
we could in principle avoid the use of our formula at this stage and instead simply 
calculate ||/||^(„+i)/ 2 (Rn) directly for the extremiser, but we wish to try to identify 
precisely the contributions of the boundary of the ball to its magnitude, and in 
any case application of the formula is fairly direct.) Leinster and Meckes m have 
recently developed an alternate approach to calculation of ||/||^(„+i)/ 2 (r„) which 
bypasses our Theorem [51 but which does not emphasise the role of the boundary. 
See the remark at the end of Section |5l 

In fact, the explicit identification of the extremiser is not an entirely straightforward 
process. Our procedures amount to giving an algorithm for its identification in any 
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odd dimension, and thus yield an algorithm for the formula for the magnitude of a 
closed ball in any odd dimension. The formulae become more complex as dimension 
increases and so we provide explicit versions only in dimensions 1, 3, 5 and 7. 

The first upshot of these calculations is that the Leinster-Willerton convex magni¬ 
tude conjecture is true for closed balls in three dimensions: 

Theorem 2. The magnitude of the closed ball of radius R > 0 inR^ is 

r3 

— + + 2R+1. 

However, in higher odd dimensions the magnitude is not a polynomial in the radius, 
thus disproving the convex magnitude conjecture in general. In dimension five we 
have: 

Theorem 3. The magnitude of the closed ball of radius i? > 0 m R® is 

R^ , 3R^+ 27R'^+ 105R^+ 216R^+ 72 
^ 24(i?-h3) ■ 


Note that when R = 0 this expression takes the value 1, and as i? —>• oo the leading 
term is i?^/5!, but the formula does not agree with the conjectured value of 


5! 




2i?2 + ^ + l. 
o 


If one expands the expression for the magnitude asymptotically for i? > 3 one 
obtains 





17R^ 21R 




as i? —oo. (Notice that some, but not all, of the coefficients here differ by a 
multiplicative factor of 9/8’ths from the conjectured values.) The role of the value 
R = —3 remains mysterious. 


Although the magnitude is not a polynomial for odd n > 5, it is the next best 
thing: 


Theorem 4. The magnitude of the elosed ball of radius R > 0, in a Euclidean 
space M" of odd dimension, is a rational function of R with rational coefficients. 


It seems possible that one might be able to take the coefficients in the polynomials 
featuring in the rational function to be nonnegative, but this has not been verified 
in general. One can show that the denominator can be taken to be a polynomial of 
degree strictly less than • 

There remain the tantalising questions of whether magnitude is a valuation in 
dimension three - and indeed whether the convex magnitude conjecture is true in 
general in dimension three - and of what happens in even dimensions, especially 
dimensions two and four. Is |i?A| a rational function of R for general compact 
convex bodies X in odd dimensions? If so, do the coefficients of the polynomials 
have a geometric significance? (See Section [TO] for a brief discussion of cuboids and 
ellipsoids in three-dimensional space.) 
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1.4. Structure of the paper. The paper falls into three main parts. 

In the hrst part, in Section[5J we begin by giving the asymptotic results of Theorem 
HJ based on little more than elementary Fourier analysis. 

In the second part we develop the general PDF theory for our problem. In Section 
Owe set out the basic facts about Sobolev spaces which we shall need. In Section 
|4] we formulate the variational problem for whose solution we develop a formula in 
Section [SJ 

In the third part we focus on the spherically symmetric situation. In Section |6] we 
give some combinatorial preliminaries. In Section [7] we find the general solution of 
the PDF problem from Section |4] in the spherically symmetric case. In Section [8] 
we fit the boundary conditions and develop an algorithm to identify our explicit 
solution and its magnitude. The discussion here leads to the proof of Theorem [4] 
in Subsection 18.41 In Section |9] we implement the algorithm in dimensions I, 3, 5 
and 7, leading to Theorems [5] and O 

Finally, in Section [TU] we make some concluding remarks. 


1.5. Acknowledgements. The authors would like to thank Tom Leinster for in¬ 
troducing them to the fascinating and beautiful subject of magnitude, and for pa¬ 
tiently and carefully explaining its category-theoretic origins in a series of seminars 
and other discussions at the University of Fdinburgh. They would also like to thank 
Mark Meckes for a number of illuminating email exchanges and comments. Finally, 
they would like to thank the referees whose careful reading of the manuscript has 
led to several clarifications and presentational improvements. 


2. Asymptotic results 


We prove the asymptotic statements of Theorem [T] which are based on elementary 
Fourier analysis. 

Let us first consider the asymptotic behaviour as R ^ oo. If / G 

satisfies / = 1 on a compact set A, it must manifestly satisfy ||/||^(„+i)/ 2 (j{n) > 


li 2 (Rn) > Vol(A), so that 


(7) 


1^1 > 


Vol(A) 


nlUJn 


as has been observed by Leinster m and Meckes m- But also 


= i?” 




so that, by the monotone (or dominated) convergence theorem. 


(8) i?-”||/(i?-'-)llff(.+i)/2(K.) 

as i? —?> oo. So we have 


l/(0PdC=/ |/(x)pdx 




Vol(d^) < ^ ■)llg(n + l)/2(Rn) ^ /gn \f{x)\‘^dx 

nlujn ~ ~ 


n\uJn 


nluJn 
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using Q, ([2]) and ([U successively. Now we can find / G with / = 1 

on X such that \\f \\2 is as close as we like to Vol(X). Indeed, with $:!&”—>• 
[OjOo) a smooth function of compact support in {|a:| < 1} with / $ = 1, $r.(a:) := 
and Xr := {x G R" : d{x,X) < r} we have that fr ■= $r * is 
a nonnegative smooth function of compact support which satisfies fr{x) = 1 for 
X G X, and /r(x) 0 as r —5> 0 for x ^ X. So by the dominated convergence 

theorem, J ^ Vol(X) as r —0. Therefore 


as i? —?> oo. 


R-^\RX\ 


Vol(X) 


Now let us consider what happens as i? ^ 0. Note that we may assume that n is 
odd since as we remarked above in Section [U magnitude is intrinsically defined. 

If / G and /(O) = I, then 

1 = |/(0)| < sup|/(x)| < ||7 ||li(B") 


so that ^ 

||/||^(n+l)/ 2 (Rn) > (^J _^47r2|^|2)(n+l)/2) 

with equality if and only if /(^) is the scalar multiple of (I + 47r^|f with 

/(O) = I, or, equivalently, /(x) = By (jS]) and Fourier inversion. 


so that 


d? 


(1 +47r2|^|2)(’^+l)/2 


I 

nlljJn 


e-H(Ode=^e-l°l 

nlOJn 


I 

n!a;„ 


(9) ll/ll J^(n + l)/2(Rn) > nluJri 

(with equality if and only if /(x) = /o(x) := This implies that |X| > I so 

long as X ^ 

Now, for 0 < i? < I, let /r : K” —^ [0,1] be a smooth function which satisfies 

1, |x| < R 

gflg-|x|, l^-l > ^1/2 

and, for i? < |x| < and |a| > I, 

d 


fnix) = 


dx 


fnix] 




We then have, (see O), 


ll/flil ^(n + l)/ 2 ( 


(n+l)/2 ^ n+1 

= E 


lD-’fi{(x)l^dx 


(n+l)/2 ^ n+1 

= E 


\ j J J\x\<R^/^ 


(n+l)/2 

\DifH{x)\^dx+e^^ Y, ' " 


\ J / J\x\>R^/^ 


l-DVo(a;)pdx 
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= / + //. 

We can estimate / by the term corresponding to j = {n + l)/2, that is, 

/ < j^n/2 _ 

and by the dominated convergence theorem, II tends to ||/o||^(„+i)/ 2 (Rn) as i? —?> 0. 
Therefore 

||/R||//(n + l)/2 ||/o||//(n + l)/2 = TlluJn 

as R —>■ 0. This shows that |i?(0, i?)| — )• 1 as i? — J- 0. 

Since for any nonempty compact set X we have (after suitable translation) {0} C 
RX C B{0,RM) for some M > 0, we immediately deduce that \RX\ —>■ 1 as 
i? —?> 0. (We thank Mark Meckes for pointing out this last implication to us.) 


3. Preliminaries on Sobolev spaces 


Let fl be an arbitrary open set in K”. We need to consider some variants of the 
Sobolev spaces consisting of complex-valued functions whose weak deriva¬ 

tives of order up to and including m G N belong to L^(n). We shall assume the 
standard properties of iL™(K") and iL’"(Ll) without specific mention. See [5] or 
[18] for more details. 


Firstly, the space iL™(n) is the completion of (fl), the class of smooth functions 
of compact support inside ft, under the inner product 


( 10 ) 





D^f-D^g 


where f 
throughout 


for j even and f = for j odd. Here and 


V/:= 



dxl 


The inner product (•, ■)h^ gives rise to the norm || • \\h^ and when = R” we can 
easily see using the Fourier transform that 




We shall be working from now on with real-valued functions, and so we shall sup¬ 
press the complex conjugate occurring in (uni). 

Let K C R" be compact and convex with nonempty interior. If / G iL'"(R"') and 
g G iL""(intiF) we say that / = g on AT if ||/-5||^m(i„tR) = 0. If / = g on AT in this 
sense, then clearly f = g almost everywhere on int K; conversely if / G iF'"(R"), 
g G iL"‘(intiF) and f = g almost everywhere on intK, then the weak derivatives 
of / and g of order up to and including m on intK coincide, and so f = g on K. 
If m > n/2 the Sobolev embedding theorem implies that the functions in //'"(R") 
are continuous and hence make sense pointwise on K\ so the statement f = g on 
K can be interpreted pointwise in this case. 

With this in mind, for our second variant we define 

iLo"*(R" \ K) ■= {/ G Fr'"(R'*) : / = 0 on K}. 
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rm /TD>^ 


), and there are natural isometric em- 


This is clearly a closed subspace of 
beddings 

\K)^ \K)^ 

Under our hypothesis that the interior of K is nonempty, the first embedding is 
surjective: 

Lemma 1. If K has nonempty interior then for any f € 77"* (R") such that / = 0 
on K, and any e > 0, there is a (j) G \ K) such that \\f — < e. 

Proof. See for example Theorem 5.29 of [2]. □ 


We need some more lemmas. 

Lemma 2. If f is in i7™(R”), K is any compact subset o/R" and e > 0 we may 
choose a smooth cut-off function if which is identically 1 on K and is such that if 
f = iff, then 11 / - /llffm(jjn) < e. 

Proof. If $ is a standard normalised bump function which is identically one on the 
unit ball, and zero off the ball of radius 2 we have 


11/(1 

and for / in 77™ (R") we have 


7=0 


'\x\>R 


|77VI' 


E 

7=0 


'\x\>R 


l^^/P —S' 0 as 7? —7 cx). 


□ 


Lemma 3. Suppose that f G 77™(R") and that / = 1 on K where K has nonempty 
interior. Suppose that if is a smooth cut-ojf function which is identically 1 on a 
neighbourhood of K. Let fix) = (/(x) — l)^(x). Then f G 77™(R” \ K). 

Proof Clearly / belongs to 77™(R”) with ||/||ffm(Rn) < C'(||/||//m(Rn) + 1) where 
C depends on if and K. Moreover / = 0 on 77. So by Lemma [2 / belongs to 
77™ (R" \ 77) with the same control. □ 

For convenience we also recall the facts about traces we shall need from [2] or [ 5 . 
(See also [T].) In the lemma which follows, denotes any differential operator 
of order j with constant coefficients. (We shall apply it in the special case that 

T)7 = foj. 

j even and = {v ■ for j odd, where v is some unit 

vector.) 

Lemma 4. Suppose f G 77™(R") and that 1 < j < m — 1. Then 

(i) If S is a piecewise smooth compact hypersurface, then T>b f belongs to 77™“-!“^ 
and 

llPVIsllff'.-.-U.(S) <C'||/||ff™(Rn); 

moreover the map f 1—7 'D^ f\s varies continuously with small changes in S. 

(ii) 7/77 is a compact convex set in R" with nonempty interior and f G 77™(R"\ 
77), then f\aK = 0. 
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Proof, (i) Since / € we have G (R") (with norm control) 

and so by the classical trace inequality its restriction to the hypersurface S is in 
provided m — j > 1, again with norm control. To see the continuity, 
suppose S and S' are nearby hypersurfaces, and take (f G and e > 0 such 

that 11/ - ())||//r„(Rr.) < e. Then 

\\V=f\s-V^f\s'\\H^-^-u^ 

The first and third terms are dominated by e by the first part and the second term 
goes to zero as S' ^ S since (j) is smooth. 

(ii) This is a consequence of the classical characterisation of the traces of functions 
in the Sobolev spaces \ K). Alternatively, it follows from the previous 

part since for smooth hypersurfaces contained in the interior of K we will have 
= 0 for all j > 1. 

□ 


4. The variational problem for the extremal energy 

In this section we consider the variational problem 

inf{||/||j^™(Kn) : / G f = g on K} 

where 5 is a prescribed member of i/™(]R”), to G N and AT is a compact convex 
set in R" with nonempty interior. We call this problem the m-extremal energy 
problem for K, or simply the extremal energy problem. 

Proposition 1. The extremal energy problem 

inf{||/||H™(Rn) : /Gif™(R"), / = 5 on K} 

has a unique solution. 

Proof. This is just the elementary fact that in the Hilbert space there is 

a unique closest point to zero in the nonempty closed convex (in fact affine) set 
{||/I|h-(R") : / e Af™(M"), / = 5 on i^}. □ 

When if is a Euclidean ball we can see immediately that the unique solution to 
the extremal energy problem is radial since averaging any solution over rotations 
yields another solution which is radial. 

By standard arguments (see for example m for the case to = (n + l)/ 2 ), the 
Euler-Lagrange equation for the extremal energy problem is 

(/ - A)"*/ = 0 on R" \ a: 

in the weak sense (testing against functions in C)(“(R" \ K)) and so we are led to 
study the problem 

(/ — A)"*/ = 0 on R" \ A' in the weak sense 
f = g on K. 

We immediately have existence of solutions in i7'"(R") to this problem from Propo¬ 
sition [3 and we also have uniqueness: 
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Proposition 2. Suppose that g G Then there is a unique solution to the 

problem 


( 11 ) 


(/ — A)*”/ = 0 weakly on R." \ K 
f = g on K 


with / e i?™(R”). 

Proof. The space of solutions to m is /o + At where /o is any particular solution 
and M is the orthogonal complement of i?^(R” \ K) in t?^(R” \ K). But by 
Lemma [TJ A4 = {0}, and so the solution is unique. 


□ 


So the unique solution to the extremal energy problem of Proposition [T] is also the 
unique solution to (fTTll in iJ'"(R"). 

We shall need in the next section to appeal to the theory of elliptic regularity. 
Briefly, for an elliptic operator such as (/ — A)"*, weak solutions / to (/ — A)™/ = 0 
on an open set are actually smooth, and the equation (/ — A)™/ = 0 holds in the 
classical sense. See [2] or [S]. 


5. A FORMULA FOR THE m-EXTREMAL ENERGY 


In this section we develop the promised formula for the m-extremal energy. See 
Theorem [5] below. 

We first consider arbitrary solutions h G iL"*(R") of the problem 
( {I- A)™ h = 0 weakly on R" \ A 

1 h = 1 on K 


where m > 1 is an arbitrary integer and A is a compact convex set with nonempty 
interior. Consider a regularised distance function d for A, which is defined on R" 
and satisfies d{x) ^ dist(a;. A). Let A^ = {x G R" : d{x) < r} and let v denote 
the unit normal pointing out of A^. (Note that dK^ is smooth and so Lemma 0] 
will be applicable.) 

For g G A'"(R"') with compact support we have 





( 12 ) 



by the dominated convergence theorem. 

We shall study the terms 9 ■ D^h by integrating by parts. In so doing, 

we shall systematically use Green’s formulae 


and 
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(where we are trying to increase the differentiability of ip, which will be a function 
of h, with the eventual aim of using the equation satisfied by h, and decrease that of 
(p, which will be a function of g). In these formulae the region will be {d{x) > r} 
and u is the unit normal pointing out of Kr. The smoothness of h on coming 

from elliptic regularity, and the compact support of g will ensure that each use of 
these formulae is valid. In particular the compact support of g means that there 
are no boundary terms at infinity to consider. 

To carry out these calculations it is convenient to define f := j for j even 
and j ■= for j odd. (Note the distinction here when j is odd 

between I) and D. While a risking a possible typographical confusion, at a first 
reading one should simply think of and as denoting suitable derivatives of 
order j. Note also that these boundary operators occur in the PDE literature in 
problems involving hybrid Dirichlet-Navier boundary conditions. See for example 
[B], p.33. It is in the lemma below that expressions such as Au (rather than 
appear more naturally on the boundary. 

Lemma 5. Let K C R" be compact and convex. For 1 < j < m and g S H^{W^) 
with compact support we have 


i-i 


[ D^g-D^h= {-ly i [ gA^h + ^{-iy [ V'^gV^^-^-^hds] . 

J{d{x)>r} \a {d{x)>r} JdKr j 

Proof. Let us first consider the case j = 1. Then 
\7g ■ Vh = — 




[ gAh- j g^dS 

J{d(x)>r} JdKr 

( f gAh+ [ gVhds] 

\y {d{x)>r} JdKr J 


as required. 

Now, for j = 2 we have 


' {d{x)>r} 


AgAh = - [ Vg- VA/i - [ ^AhdS 

J{d{x)>r} JdKr 

= ( gA^h+ f g^dS- f ^AhdS 

J{d{x)>r} JdKr- OV OV 

= [ gA^h+ [ gV^hdS- [ VgV^hdS 

J {d{x)>r} JdKr JdKr 


as required, noting that the term ^ makes sense as a member of of the 

hypersurface dKr by the trace inequality, Lemma |H (For this case to occur we 
must have m > 2.) 
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Similarly, for j = 3, we have 
I \/Ag-VAh 

J {d{x)>r'\ 

= - f AgA'^h- f Ag^^dS 
J{d{x)>r} JdKr- 


' {d{x)>r} 


Vg ■ VA^/i 


' {d(x)>r} 


gA^h- 


f dA'^h 


A'^hdS- 
dS+ f 


f A dAh 

IdKr 


dg 


A'^hdS- 


' {d{x)>r} 


gA^h+ / gV^hdS- 


IdKr 


JdKr 

f VgV^hdS^ 

IdKr 


IdKr 


A 


V^gV^hdS 


IdKr 


noting that the terms ^ and Ag make sense as members of and 

respectively of the hypersurface dKj. by the trace inequality, Lemma SI (For this 
case to occur we must have m > 3.) 

Continuing, we see that the general case takes the form 


= (- 1 )^ 


{d{x)>r} 


gA^h 


' {d(a:;)>r} 

dAi-^h 
9- 


D^g ■ D^h 


JdK, 


IdKr 

/ ^9 —^-dS'- 

IdKr 


ou 

9^9 ^,-2 


, „ -A^-2/id5 + ...| 

JdKr '* 

where the last term inside the curly brackets is 

_ / o;a - 9^j/2f^^g 


IdKr 


dv 


when j is even and 


when j is odd. That is, 


JdKr 


{d(x)>r} 


aA(j-d/2/i 

dv 

D^h-D^g 


dS” 


= (- 1 ) 4 / gA^h+ [ gV^^^-^^hdS- [ VgV^'^^-^'^hdS 

^ J{d{x)>r} JdK-r JdKr 

+ [ V^gV''^^-^'>hdS - [ V^gV^^^-'^'^hdS 

JdKr JdKr 

+ + V^^-^^gV^hdSy 

Note that the terms 'Dg.'D^g,... ,'D^^~^^g make sense as members of 
jjm-5/2^ ..., respectively of the hypersurface dKr by the trace inequal¬ 

ity, Lemma SI (For this case to occur we must have m > j.) This establishes the 
lemma. □ 
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If we take the identities of Lemma [SJ multiply by (™) and sum from j = 1 to to we 
get 


E 


\J / {d{x)>ry 


D^g-D^h 




' {d{x)>r} 


U'=l 


i=i 






k=0 


dKr 


Now add 9^ (corresponding to the the term j = 0) to both sides to obtain 


El'") / D^h-D^g 

j=0 ^ •^{d(x)>r} 

[ 9(f:h)i-iyA^h]+ Y: 'D'^g'D^^-'^-^hdS 

J{dix)>r} \^VJ/ I 0<^j<m \jJJdK^ 


{d{x)>r} 


g(/-A)™h+ V [ V^gV'^^-'^-^hdS 

0<k<j<m \g / JdKr 

y f V^gV^^-^-^hdS 

\j J JdKr 


0<k<j<m 

since h satisfies the equation (/ — A)'^h = 0 in the classical sense on {d{x) > 0} by 
elliptic regularity. 

Combining this with (1121) . we obtain the representation 

.y,. (7) L 

valid for any g € iL’”(K.") with compact support. 


We would like to extend this formula to be valid for all g € iL”^(R”). Observe that 
all the terms on the right hand side make sense for such g - as was noted during 
the proof of Lemma [5] In particular there is no problem with the existence of the 
limit as r 4- 0 for such g - simply multiply it by a smooth cut-off which is identically 
1 on {d(a;) < 2} to obtain g of compact support for which the corresponding limit 
exists and the value of which is left unaltered by the multiplication. 


So the linear functional A given by 

A 5 :=( 5 ,h)ff^(R„)-/ g-lim( ^ (-iy+>^Mf V>^gV^^-^-^hdS 

Jk yo<k<j<m \J ) j{d{x)=r} 

is well-defined on //’"(R") and is identically zero on the dense subspace consisting 
of functions of compact support. If g is in 7J"*(R"') and e > 0 we may choose by 
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Lemma [2] a cut-off function ijj which is identically 1 on {d(a;) <2} and such that if 
g = V'5, then - g\\H^(R^) < £• Hence 

^g = ^g- ■^g = Hg -g) = {g-g, h)^m(Rn) 

so that 

lA^I < ||g -5||ffm(R~)||h||/i-m(Rn) < £ ||11//m (Rn ) . 

Thus Ag = 0 for all g € and we have: 


Proposition 3. For m € N, K any compact convex set in R" and any solution 
h € iJ™(R"’) to the problem (I — A)™/i = 0 weakly on R" \ K, we have the formula 


(13) {g,h)H^(K^) 



V'^gV'^^-^-^hdS, 


valid for every g £ 


It is intuitively reasonable that for functions g £ H^fRF \ K) there should be no 
boundary terms in formula p3p and indeed for such g it is natural that g = 0 
and also that {g,h)H^ =0 using the equation. This we establish next. (Note that 
in the following lemma we had better use the equation satisfied by h as otherwise 
there is no reason to believe the terms involving higher derivatives of h exist.) 


Lemma 6. For g £ iJ™(R"’ \ K) and h £ iL™(R"’) any solution to (/ — A)™/i = 0 
weakly on R" \K we have 


lim 

r4,0 


\jJJdKr 


0<k<j<m 


V^gV'^^-'^-^hdS = 0 . 


On the other hand, if g £ iL"*(R"), if g = 0 on K (i.e. g = 0 pointwise on K if 
m > n/2, and is zero on int K when m < n/2) and if 


lim 

r .|0 


E 


I m 


IdKr 


V^gV'^^-'^-^hdS = Q 


for every solution h £ Ff'^fRF) to (I — A)™/i = 0 weakly on R" \ K which also 
satisfies h = 0 on K, then g £ iL™(R" \ K). 


Proof. Consider the first assertion. By our formula (1131) we have 

hm E (-1)'+"W / V'^gV^^-'^-^hdS={g,h)H^ 

\3jJdKr 


0<k<j<m 

So it is enough to show that 


'K 


Ag := {g,h)Hm(R,^) - g 

JK 


satisfies 

A5 = 0 

for all g £ H^(R" \ K). Let e > 0 and g £ C^(R" \ K) be such that ||(; — 
5 ||jp™(R"\a:) < £• Clearly we have A^ = 0 since h satisfies (/ — A)™/i = 0 weakly 
on R" \ iL and if n R" \ if = 0, and so 


- / {g-g)h, 

JK 


^g = ^g- ^g = Hg -g) = {g-g, ^)ff-(R") 
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thus 

IA 5 I < II 5 - + \\g - g \\2 ||/l ||2 < £ 11 /l 11 //m (Rn ) . 

Hence Ag = 0 for all g G \ K). 

For the second assertion, if the limit in question is zero, we have by the represen¬ 
tation formula (USD that {g,h)H^ = 0 for all h G which are solutions to 

(/-A)™/i = 0 weakly on R”\iF which also satisfy h = 0onK. So g G 
and g is perpendicular to the space A4 introduced above in Proposition [21 By def¬ 
inition of A4 this means that g G \ A"). 

□ 


We have not yet used the requirement that AT have nonempty interior nor used any 
boundary conditions satisfied by h (except to obtain the formula (HID)- In order 
to deduce a useful formula for ||Ii||^m(Rn), we shall have to do so, and we shall use 
Lemma [3] freely from now on. 


So, let H be any function in iL'"(]R”) such that H = 1 on K. Then, by Lemma 
IS] we have that g{x) = {H{x) — l)'ip{x) belongs to \ K) (where i/i is any 

smooth function of compact support which is identically 1 on {d{x) < 2 }), and so 
we can apply Lemma [ 6 ] to obtain 


lim (- 1 )^+''/'™') f V^HV^^-^-^hdS 

0<k<j<m ^ JdKr 

-liin Y = 0. 

0<k<j<m JdKr 


But for r sufficiently close to zero we have 


so that 


(14) 


\ J / JdKr 






= 11 


'dKr 


lim Y f V^HV^^-’^-^hdS 

\J / d dKr 


0<k<j<m 
m 

= lim V(-1)J 

rJO ^ 

1=1 


V^^-^hdS. 


J J JdKr 


Next, we take g = H in Proposition |S] and apply (HTl) to obtain: 

Proposition 4. Suppose AT has nonempty interior, h is any solution in iL’"(IR."\iL) 
to (/ — A)™/i = 0 weakly on R" \ K, and H is any function in iL’"(IR.") such that 
H = 1 on K. Then we have 

{H,h)Hm^R.)=YoliK)+limYi-iyh) f V^^-^hdS. 

\jJJ 0 Kr 
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Finally, we claim that the terms with j < mj2 give a contribution of zero in this 
formula. 


Indeed, if we have h G and h = 1 on K, Lemma [3] tells us that for 

suitable smooth ip of compact support, identically 1 on a neighbourhood of K we 
have g = [h — Vjip G \ K). So by Lemma |4l V^g is the zero member of 

for 0 < j < m — 1, and continuity of traces shows that we therefore 


have 



V^gdS —5> 0 as r J, 0 


for 0 < i < m — 1. 


In particular, since V^h = V^g + near dK we have 

I V^hdS —>■ 0 as r 4- 0 

JdKr 

for 1 < i < m — 1. So 


f V^^-^hdS ^ 0 as r ; 0 
JdKr 


for 1 < — 1 < m — 1, that is for 1 < j < m/2. 


Hence 

=Vol(iL)+lim V i-iyh) [ V^^-^hdS, 

and we have proved: 


Theorem 5. Suppose that K is a compact convex set in R” which has nonempty 
interior. Let H be any function in H^(W^) such that H = 1 on K. Then the 
unique solution h G to the problem 

(/ - A)^h = 0 weakly onR^\K 

h = 1 on K 

satisfies 

{H, = Vol(iL) + V (-ly W lim [ -^A^-^hdS. 

In particular, taking H = h, 

(15) ||/i|ll,™(r)=Vol(iL)+ ^ (-l)^Mhm/' ^A^-^hdS. 

/o ^ __ \ 3 / J dKr 

ml2<j<m ^ 


If the boundary of K is sufficiently regular to allow the invocation of boundary 
regularity for elliptic equations we may realise this formula more succinctly as 

The main virtue of Theorem [5] is that it explicitly demonstrates the dependence of 
the TO-extremal energy upon the volume of K and integrals over its boundary. It is 
analogous to the representation of the classical capacity inf{/ |Vup :M=loniL} 
as for u satisfying Au = 0 off Lf, u = 1 on dK and u ^ 0 at oo. The 
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fact that = (/i, for any H G which is identically 1 on 

K when m = (n + l)/2 is more general and due to Meckes [T7], Proposition 4.2. 
One may wonder whether the alternating sum in Theorem [5] reflects an underlying 
“Euler characteristic”, and whether the individual terms appearing might have 
some separate meaning. It is also of interest to determine whether or not the 
solutions of the boundary value problems are always positive. See for example [6] 
for a discussion of such matters. 

As mentioned in the Introduction, Leinster and Meckes have an alternative ap¬ 
proach to calculating magnitudes. Indeed, Theorem 4.16 of [13], specialised to the 
case of Euclidean spaces, shows that up to a constant, the magnitude of a compact 
subset of R." is the integral of its potential function, provided the potential function 
is in L^(R"). In the context of Theorem^and formula dTS]) this translates simply 
as 

= f h. 

JR” 

While this formula arguably offers some comptational advantages over formula ([TK]) , 
it does not bring out the important role that the boundary of K plays. In this regard 
see also the last remark in Section [101 

6. Combinatorial preliminaries 

For our discussion of solutions to certain ordinary differential equations we shall 
need to consider some special polynomials which have a combinatorial flavour. 

We define goit) := 1 and, given gj, we define gj+i by 

9j+i{t) :=t^gj{t)+tgjit). 

Lemma 7. For j > 1 we have 

9j(t) 

k 

where = 0 for k < j and k > 2j, = 1 for all j, and each cf for j < k < 2j — 1 

is a nonnegative integer. 

The easy proof is left to the reader. The defining formula for gj leads us immediately 
to a recurrence relation for the coefficients c^: 

Lemma 8. For j + l<k<2j — 1 we have 

(16) <4+1 = (fc - l)c4_i -b cf 

and 

(17) 4l'+i = (2j-l)4,-i- 

Moreover for J T 1 < fc < 2j — 1 we have 

2j<44i = (fc - 1)^4-! + 2fc'4- 
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Proof. We have by definition 

V fc=i / V / 

= E + E = E + E 

k=j k=j k=j+l k=j 

fc=j+i 


On the other hand, 


2i+l 2j 

<7j+iW = E 

fc=i+l k=j 

and then (HU) and (II3 follow by comparing the two expressions. The final identity 
follows readily from the expression (ITSl) below. 

□ 


Pictorially, identity (ITSl) can be represented as a Pascal triangle as follows, (with 
the first row representing c\, the second c^, c|, etc.): 

1 

1 1 
1 3 3 
1 6 15 15 
1 10 45 105 105 
1 15 105 420 945 945 


An explicit formula for for j < k < 2 j — 1 is given by 

^ _{k-l){k-2)---{2j-k) 

2^-^{k-j)\ 

and this follows immediately from equations (HH) and (ini) together with the fact 
that cj = 1 for all j. 

The numbers c(, are reminiscent of Stirling numbers and are sometimes called Bessel 
numbers of the first kind; the gj are related to Grosswald’s polynomials. See [8] 
and [21j . The gj can also be expressed in terms of generalised Laguerre polynomials 
and hypergeometric functions; we thank Jim Wright, Chris Smyth and Adri Olde 
Daalhuis for pointing these connections out to us. 
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7. Explicit radial solutions to (I - A)™/i = 0 

When if is a Euclidean ball, the unique solution to problem (HU with 5 = 1 is 
necessarily radial (since an average over rotates of a solution is also a solution). So 
in this section we seek radial solutions h G iJ"*)®”) to 

(/-A)™/i = 0on{|a;| > R} 

where R > 0. Then we shall impose the extra (boundary) condition /i = 1 on {|a;| < 
R}. We are particularly interested in the case where n is odd and m = (n + l)/2. 
We shall systematically abuse notation and not distinguish between the function 
h{\ ■ I) defined on {|a;| > R} in R" and the function h{-) defined on {R, oo). 

Recall that the action of the Laplacian in R" on radial functions is given by 

rf _ 1 

A/(r) = /"(r) + - - f{r) 

r 

and so we define, for an integer u > 0, 

A./ := r + —f 

r 

and consider solutions of 

(19) {I-A,rh = 0 

on (i?, oo) where R> 0, with h{\ ■ |) S L^(K"). (We hope that the the use of u both 
as a proxy for (n — l)/2 in this section and as a normal direction to K in previous 
sections will not cause confusion.) If n is odd and v = (n — l)/2 this amounts to 
finding radial solutions in of 

(/- A)”"/! = 0 on |a;| > R. 


Let us set up some notation. We define, for A: > 0, 


fk{r) 


We define = e ^ and, for j >1 


e 


—r 


r 


k 


2j-l 

k=j 

where the gj are defined in Section [6] 

From Pascal’s triangle for the coefficients c^, (Lemma [5]), we can write down the 
first few ipj explicitly: 


^o = fo = e-\ 

'ipi = fi = Ir, 

V'2 = /2 + /s = e“’'(r“2 + r“^), 

i’3 = f3 + 3/4 + 3/5 = + 3r“'‘ + 3r“®), 

i/’4 = /4 + 6/5 + 15/6 + 15/7 = e“’'(r“^ + 6r“® + 15r“® + 15r“’^), 

^5 = /s + lO/g + 45/7 + 105/8 + 105/9 = e"’'(r“® + lOr"® + + 1057"® + lOSr"®) 
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Lemma 9. For k >0 and v >0 we have 

- I)fk = 2{k - v)fk+i +k{k+l- 2 v)fk+ 2 - 

In particular, 

(A,-/)/o = -2u/i. 

Proof. Let us calculate 


so that 


fk = = -fk - fc/fe+i 


fk = -fk - kfUi 

= -{-fk - kfk+i) - k{-fk+i -{k + l)fk+2) 
= fk + 2kfk+i + k{k + l)/fe+2- 


Therefore 


(A, - I)fk =fk + 2kfk+i + k{k + l)fk+2 - 2u f ^ ) - fk 

\ r r ' 

= 2kfk+i + k{k + l)fk +2 - 2v {fk+i + kfk+ 2 ) 

= 2{k - i')fk+i + k{k + l- 2i2)fk+2- 


Proposition 5. For j > 0 and v >Q we have 

(A^ - I)ipj = 2{j - 

In particular, 

(A, - I)^, = 0. 

Proof. The case j = 0 is a special case of Lemma [9l For j > 1 we have 

2i-l 2j-l 

(A. - I)^|;j = (A. - I) ^ cifk = 4 (A. - I)fk 

k=j k=j 

2j-l 

= X! 4 [“2(u - k)fk+i +k{k + l- 2iy)fk+2] 

k=j 

by Lemma [9l Rearranging, this equals 
2j-l 2j 

2(* - v^kfk+i + ^ {k- l)(fc - 2u)<4_i//c+i 

k=j k=j+l 

= Hi - v)c’jfj+i + (2j - l)(2j - 2 u)<4^_i/2j+i + 

2j-l 

+ {2(/c-u)<4 + (fc-l)(fc-2i/)<4_J/fc+i 

k=j+l 

= 2(j - v){fj+i + {2j - l)4^_j2j+i}+ 

2i-i 

+ {2(fc-i^)<4 + (^-l)(^-2i^)4-i}/fc+i- 

k=j+l 


□ 
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Now since 

2i+l 2j 

2(j - v)'tljj+i = 2{j -v) = 2(j - u) ^ 4+i/fc+i(^) 

k=j+l k=j 

we see that the coefficients of /j+i agree (as c^- = 1 for all j), and it is therefore 
enough to see that 

4^+1 = (2j - l)4j-i 

and, for j + 1 < fc < 2j — 1 and arbitrary v 

2(j - ^)ci+i = 2(fc - + (fc - l)(fc - 2u)c4_j. 

This latter splits as 

4+1 = 4 + (fc - 1)4-1 

(coefficient of v) and 

2j44i = 2^4 + (fc - l)fc4-i 

(u = 0). All of these are true identities by Lemma [S] □ 

Corollary 1. For j >Q we have 

( 20 ) ip'jir) =-ril^j+iir). 

For m > 1, V > 0, and either 0 < j < v — m or j > v we have 

(21) {A^ - /)”"V'j = 2™(j + m - 1 - u)(j + m - 2 - u) • • • (j - 
while for v — m < j < v we have 

(22) (A, - irfjj = 0. 

Finally, for m > 1, v > 0 and ip = El=o 

y—m m —1 / \ 

(23) A^fip = ^ a,2-(j + m - 1 - u) • • • (j - ^ (-1)—'= 

Proof. Write the first identity of Proposition 0 as 

4' + “ 4' = 2(j - T^)'ipj+i, 

divide by u, and let u —>■ oo to obtain (uni)- 

Let j > 0 and m > 1. By repeated application of Proposition [5] we obtain (|^ and 
([22]). 

By the binomial theorem, 

m ^ V 

(A. - in = 

so that for j < u — m, by (EU, 

m-l , , 

(24) Affip, = 2™(j + m - 1 - u)... (j - v)n - E (-I)™”'' T 

k^O ^ ^ 
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while for u — TO < j < u, by ([2^ , 

m-l , ^ 

(25) = 

k=0 ^ ^ 

Finally (I23|) follows from (1241) and (|25|) by taking linear combinations. □ 

We shall use Corollary [T] to systematically turn the calculus into algebra when 
fitting the appropriate boundary conditions to the problem, and when calculating 
the right hand side of formula (ITSi) in Theorem [5] in the radial case. Note that (|23)) 
tells us that for all to > 1, if i/i = X)j=o then A^V' can be expressed in terms 

of 

Returning now to the ordinary differential equation m, we see that (EH) from 
Corollary [T] tells us that the general solution to 

{I-A,rh = 0 

on {r > 0} with h{\ ■ |) € L^(R") at infinity (for suitable n) is given by a linear com¬ 
bination of ■ ■ ■ ,4’v-m+i}- (Routine calculations similar to those carried 

out at the beginning of this section reveal that the other m linearly independent 
solutions will feature and so will not have the requisite decay at infinity to 
belong to any L^(R").) 

In particular, taking u = (n — l)/2 for n odd, the general radial solution to 

(/ - A)^h = 0 on R" \ {0} 
with h € L^(M") at infinity is 

(1^1) T * *' T (l^^l) 

where the a’s are arbitrary scalars. 

Specialising further, taking to = (n -I- l)/2, we obtain: 

Proposition 6. Suppose n is odd. Then the general radial solution to 

(/ - A)(”+i)/2/i = 0 on R” \ {0} 

with h G L^(R") at infinity is 

h{x) = aofjoilxl) H-h Q;(n-i)/2'0(n-i)/2(|a;|) 

where the a’s are arbitrary sealars. 

8. The algorithm for extremal energy and magnitude 

8.1. Fitting the boundary conditions. The next task is to determine the aj’s 
which permit the solution h = ' I) of 

(/ - A)(”+ih2/i = 0 on {|a:| > R} 
h = 1 on {|a:| < R} 

given by Proposition [5] to lie in ij("+i)/2(R"). This amounts to forcing it to satisfy 
appropriate boundary conditions at |a:| = R. Now our radial solution h belongs 
to gQ^ Lemma m its derivatives of order j lie in iJ"/2-f of 

spheres (for 0 < j < (n — l)/2) and moreover vary continuously with small changes 
in the radius of the sphere. Alternatively, the function h defined on (0, oo) belongs 
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in near |a:| = R, and so by Sobolev embedding in the one-dimensional 

case, the derivatives of h of order up to and including (n — l)/2 are continuous 
at R. In any case, since h = 1 on {\x\ < i?}, this manifestly forces the (n -|- l)/2 
conditions 

h{R) = 1 
h'{R) = 0 
= 0 

(26) {^h)'{R) = 0 


(A("-i)/4/i)(ij) = 0 or (A(’"-3)/^/i)'(ii) = 0 

depending on whether n—lorn — 3isa multiple of 4. In either case, this gives 
us (n -I- l)/2 linear conditions on the (n -I- l)/2 unknowns oq, •.. ,a(n-i )/2 which 
therefore determine them. (Since we know that there is a unique solution to the 
boundary value problem we know that the linear system (HH) has a unique solution. 
Hence the system is guaranteed to be nonsingular - otherwise there would be either 
no solution or multiple solutions.) In practice, we shall need to use Corollary [T] in 
conjunction with the explicit formulae for the ipj in order to evaluate the aj. As 
above, we write v = (n — l)/2. 

We shall use (1231) from Corollary [l] directly on the equations involving A^h(R), 
and ((^ together with (1^ on the equations involving (A^7i)'(ii). 

The first equation is simply 

aoipoiR) H-1- a„^pi,{R) = I 

while subsequently we have, by (|23ll and the boundary conditions on the terms 
h, Ah, ...(A)™-i/i, 


A'"/i(i?) = 2"* ^ (j -h TO - 1 - v){j + - v)a;j%l^j+m{R) - (-1)™. 

3=0 

Since the terms (j -|- to — r — j^) are all negative and there are to of them, the 
boundary condition A'"/i(i?) = 0 reduces this to 

v—m 

2™ ^ {v- j)(v - j - j -m+l)aj%j3j+rn(.R) = 1- 

3=0 

Similarly, the boundary condition (A'"h)'(i?) = 0 reduces to 

V—m 

2”" ^{v - j)(v - j - 1) ■ ■ ■ {v- 3 - m + l)ajij)j+m+i(.R) = 0. 

3=0 


Thus we have the -I- 1 = (n -|- l)/2 equations 
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(27) 
Av{v - 

‘iv{v - 


aoipQ{R) + aitj^iiR) 
aoipi{R) + ai'ip2{R) 

2vaQ'4’i{R) + 2{v - 1)q;i^2(-R) 
2vaoil}"2{R) + 2(j^ - l)aiV’3(-R) 
^)aoi’2{R) + 4 (z/ - l){v - 2 )Q!i^ 3 (i?) 
l)ao'!/'3(^) + 4(i^ - l)(i^ - 2)Q!iV’4(i?) 


-b-- 

+ a^tp,^{R) 

= 1 

-b" 

-b a^y^+i{R) 

= 0 

-b" 

-b 2ai/_i^[/(i?) 

= 1 

-b" 

-b 2cri/_i'0i/_|_i (i?) 

= 0 

-b" 

-b 4.2.1ai,_2V'!2(I?) 

= 1 

-b" 

-b4.2.1ai,_2V'!2+i(i?) 

= 0 


to solve for the v + 1 = (n + l)/2 unknowns ao, ai,...,ay. (We know the system is 

nonsingular, but it is perhaps not so obvious that this is the case by direct inspec¬ 
tion.) Notice that the “input” here involves only the values of iIjo{R), ..., '0j/+i(i?). 

We implement the explicit solution of these equations in the first few cases corre¬ 
sponding to n = 1,3, 5 and 7 below. 


8.2. Calculating the extremal energy. Once we have determined the aj = 
ctjiR) which satisfy (l27l) . we can use Theorem [5] to calculate the extremal energy 
ll^ll^(n+i)/ 2 (Rn) for h = aoV'o -b • ■ • -b Indeed, we have, with m = (n + l)/2. 


|^(n + l)/ 2 (Rn) — Yo\{K) -b ^ (“ 1)^ 

in+l)/i<j<in+l)/2 


n+1 


lim 


j J ’'to jQXr dv 




(28) 


{n+l)/A<i<{n+l)/2 




= UJnR'" \ I + nR-^ 

{n+l)/A<j<{n+l)/2 


\ j ) 

(T)- 


(-i)^r2'VA^-ih)'(i?) 


where LOn is the volume of the unit ball in K" and ct„_i = na;„ is the surface 
area of the unit sphere S"“^ in R". We emphasise that by (A'l“^h)'(i?) we mean 
limr 4 .R(A'^“^/i)'(r) here and subsequently below. 

The calculation of {/\^~^h)'{R) is carried out using (|2^ from Corollary[T] Indeed, 
this gives us, for j > 2, 


v-j + l 

A^~^h = "Y {i + j - 2 — iy){i + j - 3 - L') • ■ ■ {i — v)ai'ilji+j-i 

i=0 

Now invoking ipj(r) = —ripj+i{r) (that is, (1^01) from Corollary [IJ , we obtain 


{A^-^hy{R) 


iy-j+1 

R2^~^ ^ {i + j-2 - iy){ij - 3 - u) • - {i - iy)ai'ilJi^j{R) 

k^O ^ 
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The boundary conditions imply that we have (A^h)'(i?) = 0 for 0 < /c < (u — 2)/2 
when V is even, and 0<fc<(u — l)/2 when v is odd. So in the /c-sum we only need 
consider terms k with k > [v — l)/2, and therefore, recalling that v = {n — l)/2, 
we have, for j > 2, 

u-j + l 

= -i?2^-i Y. (i + j - 2 - u) • • • (i - v)aii^,+j{R) 

(29) *=° 

^ (-1)^-1-'= 

{v-l)/2<k<j-2 

This allows us to effectively compute the terms (A^“^/i)'(i?) recursively for j > 2 
in terms of the values of ipo{R),..., ip,^+i{R) and ao,..., Oi,. We need this for those 
j such that (n + l)/4 < j < (n + l)/2. Once done, we can substitute into formula 
(I28p for the extremal energy, and we carry this out explicitly in the cases n = 1, 3, 
5 and 7 below. 

8.3. Summary of the algorithm. We first solve the equations (071) for the (n + 
l)/2 unknowns oq, ..., 0 (n-i)/ 2 - Then we use formula (0^ to calculate (A^“^/i)'(i?) 
for (n + l)/4 < j < (n + l)/2, which we can substitute into formula (051) . whose 
value we then compute to give the extremal energy. Finally we divide by a;„n! to 
obtain the magnitude. 




8.4. Proof of Theorem [4l The discussion of the algorithm above immediately 
gives that the magnitude of the closed ball of radius R in an odd-dimensional 
Euclidean space is a rational function of R. Indeed, modulo the multiplicative 
exponential term e~^, all of the inputs 'iIjo{R), ... ,'ij:v{R) are polynomials in R~^ 
with nonnegative integer coefficients. The solutions to our system of 

equations will be given, by Cramer’s rule, by a ratio of two determinants, each 
of which is a polynomial in tjjQ{R),... ,tjji,(R) with integer coefficients, up to a 
multiplicative exponential term e^. Inspection of formulae ([29|) and ([28)l indicates 
that the magnitude is given by i?"/n! plus an integer combination of various arips 
for 0 < r < (n — I)/2 and 0 < s < (n -|- l)/2. So the multiplicative exponential 
terms will cancel throughout in every instance, and we will be left with i?"/n! 
plus a rational function of R with integer coefficients. 

9. Implementation of the algorithm 


tpo{r) = e 

^i{r) = Y- 

V’2(i-) = e"’' + :^) 

tpair) = 6 “’’ (tt + ;!■ + ^) 

Mr) =e-^M + ^ + ^ + ^). 


Recall that we have 
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9.1. Case n = 1. Notational fix: v — 0, m = 1. 

Inputs: ipo{R), 'ipi{R). 

Outputs: ao, h'{R), extremal energy, magnitude. 

The single equation for oq is ao'ipo{R) = 1 which has solution ao = If-ipoiR)- 

Equation (l28l) requires for j = 1 only, which is simply h'{R), and by 

Corollarylll equation ([^ . h'{R) = —Rao'fi{R). 

Using ([28|) . the extremal energy is given by 

ujiR ^1 + hfR^ = ujiR (1 - R-^h'{R)) = uJiR (1 + aoMR)) ■ 

Substituting the value of oq gives 

uJiR (1 + ipi{R)/'ipoiR)) = uji{R+ 1 ). 

Dividing by llwi gives the magnitude of the closed interval [—R, i?] as 

R 1, 

in agreement with the results of m and m- 

9.2. Case n = 3. 

Notational fix: v = 1, m = 2. 

Inputs: ifoiR), ifiiR), tjj 2 {R)- 

Outputs: ao, oi, h'{R), (Ahy{R), extremal energy, magnitude. 

Our two equations are 

aoifoiR) + ai'fi{R) = 1 

aotpi^R) + aif>2{R) = 0 

which have solution 

^ _ V^2(-R) _ 

'ii)o{R-)'f2{R-) -'fi{RY 
^ -yiiR) 

'ifo{R)i’2{R) - ipiiRy 

Now 'ipo{R)'4’2{R) — 'fi{R)'^ = e~^^{R~'^ + R~^ — i?“^) = e~'^^R~^, so 

ao = e^{R + 1) 

ai = -e^R^. 

Note the appearance of the expression i?+l (as displayed in the previous subsection) 
in the formula for ao- 

Equation (l28ll requires for j = 2 only, and by formula (l29)) we have 

{AhYiR) = -R{2)\-l)aoMR) - 0 = 2RaoMR) 

Using (1^51) . the extremal energy is given by 

W3i?3 (^1 + 3i?-i(-l)2 Q {Ah)'{R)^ = ujoR' (1 + ?,R-\Ah)'{R)) 

= ujoR^ (1 + 6ao'(/’2(i?)) ■ 
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If we substitute in the explicit values for ao and ip 2 {R) we obtain 

= uj3{R^ + 6R^+ 12R + 6), 

which upon dividing by 3!w3 gives the magnitude of the closed ball of radius R in 
as 

1 R^ 

- (i?3 + 6i?2 + I2ii + 6) = — + i?2 + 2i? + 1. 

This verifies the Leinster-Willerton convex magnitude conjecture for balls in R.^. 

9.3. Case n = 5. 

Notational fix: v = 2, m = 3. 

Inputs: i’o{R ),..., ifsiR)- 

Outputs: ao, ai, a 2 , h'{R), Ah'{R), extremal energy, magnitude. 

Our three equations are 

CIoV'o(-R) + CnV'l(-R) +Q!2'02(.R) = 1 
Q!o'0i(-R) + CH'02(I?) +a2tp3{R) = 0 

4aoi/’i(i?) + 2aii/’2(i?) +0 =1 

which upon substitution of the values of ifj (i?) become 

ao + aiR~^ +a2{R~‘^ + R~^) = 

aoR-^ + aiiR-"^ + R-^) +a2(i?”^ + 31?"^ + 31?"'^) = 0 

aoAR-^ + ai2{R-‘^ + R-^) +0 = e^. 

A visit to matrixcalc.org confirms that the solution of these equations is given 

by 

a 

^ 2{R + 3) 

where 

/3o = 6 + 12i? + + i?3 

Pi = -i?^(12 + 9i? + 2i?2) 

^2 = (2 + 

Note that the system is singular when i? = — 3, and note the appearance of the 
expression R^/6 + + 2R + 1 (as displayed in the previous subsection) in the 

formula for ao- 

Equation (1281) requires for j = 2 and 3, and by formula (l29l) we have 

(case j = 2) 

1 

(AhYiR) = -R.2j2ii - ‘2)a,Y^+2{R) + 0 = 2R{2aoMR) + ^iMR)), 

2=0 
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and, (case j = 3), 

{A^hYiR) = -R2\-l){-2)aoMR) - (- 1 )' Q {^hy{R) 

= —SRaQtpa^R) + AR{2ao'4’2{R) + o^itpsiR)) 

= —8Raoip3{R) + 8Rao'ip2{R) + iRaiipsiR). 

Using (I28L the extremal energy is given by 

u:,R^ (^1 + 5i?-i{ Q) - (3) 

= ujsR^ (1 + 5R-^{3{Ahy{R) - {A‘^hy{R)}) 

= uj^RY’ (1 + 5{3[4ao'*/'2(^^) + 2aii/>3(i?)] - [-8aotp3{R) + 8aoip2{R) + 4ai^3(i?)]}) 
= uj^R^ (1 + 5{4aoV'2(-R) + 8aoil}3{R) + 2 ai'!/)3(i?)}) . 


Substituting in the explicit values for 1/12 {R) and ^p 3 (R) we see that this expression 
equals 

wg (i?® + 5e"-^i?^[(4ao(i?-2 + R-^) + 8ao(i?“^ + 8 R-^ + 3i?“®) + 2ai{R-^ + ZR-^ + ZR-^)]) 

= wg (i?® + 5e"-^[ao(4i?3 + 12i?2 + 24i? + 24) + ai{2R^ + 6 R + 6)]) . 

Now plugging in the explicit values for ap and ai, we obtain 

((6 + 12-R + 6i?^ + i?^)(4i?3 + 12R^ + 2AR + 24) - R^{12 + 9R + 2R^){2R^ + 6i? + 6)) 

2[R + 3) 

which, after some simplification and cancellation in the highest order terms, gives 
that the extremal energy is given by 

, 72 + 216i? + 216i?2 + i05i?3 + 21R^ + ZR^ 

UJ 5 R +5wg-——-. 

ri o 

Dividing by cg = 5!wg we obtain that the magnitude of the ball of radius R in R® 
is 

R^ ^ ZR^+ 27R^+ 105R^+ 216R'^+ 216R + 72 
^ 24(i? + 3) ■ 

Note that when i? = 0 we obtain the value 1, but when i? > 0 this quantity is 
strictly greater than the conjectured value of 

^R^ + -R* + -R^ + 2R^ + -R + 2i?2 + l. 

5! 9 3 3 

We can expand our formula when i? > 3 as 

^ r5 ^ r4 I ^ d 3 I 1"^ n2 I 21 „ 9 1 

_R +_R +_R +_R +_R+_+0{-) 

as i? —!> 00 , but the coefficients (other than that of i?®) do not agree either with 
those conjectured by Leinster and Willerton. 
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9.4. Case n = 7. 

Notational fix: u = 3, to = 4. 

Inputs: '4’o{R)j ■ • ■) '4’4{R)- 

Outputs: (Toj • ■ ■; Q^ 3 j h'{R ),..., (A^/i)'(i?), extremal energy, magnitude. 

For brevity of expression we suppress the multiplicative factors which should 
appear below, and which ultimately cancel as described above. We have the four 

ao-ipoiR) + cxiipi{R) + a2i>2iR) + astpsiR) = 1 
aoipi{R) + ai'ip2{R) + oi2'4>'3[R) + = 0 

< 

6ao^i(-R) + 4aiV'2(^^) + 2a2ip3iR) = 1 

. 3aopJ2{R) + 2aiip3{R) + a2ip4{R) = 0 

_ 360+1080R+1080R^+525fi^ + 135R* + 18R^+R° 

“0 360+288R+72R^+6R“ ’ 

_ 360R^+555R^+345fi'* + 105R^ + 16R‘’+R^ 

120+96R+24R2+2R3 ; 

_ 120R^ + 150R'^+66R^ + 13R’’+R^ 

“2 120+96R+24R^ +2R3 : 

_ 24R°+27R^+9R^+R° 

“3 360+288fi+72R2+6/i3 

as another visit to matrixcalc. org confirms. 

Note that the numerator in the formula for ag is 

360 + lOSOi? + 1080i?2 + 525R^ + 135i?^ + 18R^ + i?® 

whereas the magnitude of the ball of radius R in R® (calculated in the previous 
subsection) is 

i?5 ^ 3E’5 + 27i?4 + 105i?3 + 216i?2 + 216i?+72 
^ 24(i? + 3) 

360 + 1080i? + 1080i?2 + 525i?3 + 1357?^ + ISi?^ + R^ 

120(77 + 3) ’ 

with the same numerator as oq. 

From equation (1281) the extremal energy is given by 

u:rR’^ + TiOrR^ (^- Q (A^u)' (7?) + Q (A^u)' (77)) 

^32) = ^7 {r^ + 777«(-4 (A^u)' (77) + (A^u)' (77)) 

= UJ7 (r^ - 2877® (A^u)' (77) + 777® (A®m)' (77)) 

= UJ7 (77^ + 777®{(A®u)'(77) -4(A2 u)'(77)}) . 

So we need to calculate {A'^u)' (77) and (A^u)' (77). 


equations 

(30) 

with solution 

(31) 
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Using the recurrence relation ([29| with j = 3 gives 

1 

{A^hYiR) = - 2){i - 3)a,V’i+3(i?) + 0 

i=0 

= —4i?(6aoV’3(-R) + 2aiip4{R)) 

= -8R{3aoip3{R) + aiip4{R)). 

Next, using the recurrence relation (l29ll with j = 4 gives 

(A^hYiR) = 6R2^aoMR) + Q(A"h)'(i?) 

= A8RaQip4{R) + 3{—8R{3ao'ip3{R) + aiip4{R))) 
= 48i?ao'04(^) ~ 24:R{3ao'tp3{R) + ai'ip4{R)) 

= 24i?(2aot/’4(i?) — SaotpsiR) — aiip4{R)). 


360 + 1080i? + 1080i?2 + 525i?3 + i35i?4 ^6 3 3 ^ ^2 

120 + 96i? + 24i?2 + 2i?3 ^ ^5 

360 + 555i? + 345i?2 + lOSi?^ + i6i?4 + _r5 + IhR + QR^ + R^ 

120 + 96i? + 24ii2 + 2i?3 ^ ^5 

4320 + 9405i? + 88201?^ + 4545i?3 + 1380i?4 + 246i?5 + 24i?6 + RY 
(120 + 96i? + 24i?2 + 2i?3) 

(noticing the cancellation in the i?® terms), and so 

/ A 2 N', 8(4320 + 9405i? + 8820i?2 + 4545i?® + 1380i?4 + 246i?® + 24i?® + i?^) 

(A^u) (R) = -J--^-. 

'' ’ ^ ' i?4('i20 + 96i? + 24i?2 + 2i?3) 


Now 

3aotp3 + aii’4 = 


Moreover 

3 9 9 2 12 30 30 -30 - 30i? - 3R^ + 7R^ + SR'^ 

3^3-2V’4 - ■^3+Jf4+J^-Jf4-J^-J^-^ - -’ 

SO that 

ao (SV's - 2V'4) + aiV'4 

equals 

360 + 1080i? + 1080i?2 + 525i?® + 135i?4 + 18R^ + iJ® -30 - 30i? - 37?^ + jrs + 

360 + 288R + 72i?2 + 6R3 ^ ^7 

1080i?2 + 1665i?® + 1035i?4 + 315i?® + 48i?® + 3i?^ 15 + 15i? + 6R^ + R^ 

360 + 288R + 72i?2 + qrs ^ ^7 

10800 + 43200i? + 82080i?2 + 90045i?® + 613807?^ + 26685i?® + 7380i?® + 1254i?^ + 120i?® + 5E® 
“ (360 + 288R + 72i?2 + QRS) 

(noticing the cancellation in the terms), which gives that 

{A\y (R) 
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equals 

24(10800 + 43200i? + 82080i?2 + 90045i?3 + 61380i?4 + 26685i?’5 + 7380i?® + 12547?^ + 120i?® + 5R^) 

i?6 (360 + 288R + 72i?2 + qr3) 


Combining these, we see that 

R^ (A^u)' (R) - 4i7® (A^u)' (i?) 

equals 

24(10800 + 4320077 + 8208077^ + 9004577^ + 6138077^ + 2668577^ + 738077® + 125477^ + 12077® + 577®) 

360 + 28877 + 7277® + 677® 

96(432077® + 940577® + 882077^ + 454577® + 138077® + 24677^ + 2477® + 77®) 

360 + 28877 + 7277® + 677® 

_ 43200 + 17280077 + 25920077® + 20970077® + 10440077^ + 3402077® + 744077® + 108077^ + 9677® + 477® 
“ 60 + 4877+ 12772 + 77® 


Finally, we insert this into (15^ and divide by 7!w7 to see that the magnitude of the 
closed ball of radius 77 in is given by 

77^ , + f + 14577^ + i^77® + 36077® + 24077 + 60 

TT 77® + 1277® + 4877 + 60 


10. Concluding remarks 


1. The three-dimensional convex conjecture predicts that the magnitude of the 
convex body K should be 

Vo{K) ^ Vi{K) ^ V^jK) ^ VsjK) 

0! ujQ 1! uji 2! UJ 2 3! cus 

where Vj (K) is the j’th intrinsic volume and wq = 1, wi = 2, a ;2 = tt and wa = 47r/3. 
Now VoiK) = 1, +2(7^) = Surf {dK)/2, VsiK) = Yo\{K) and Vi{K) = P 2 {K)Itt 
where 7^2(7F) is the coefficient of t® in the degree 3 polynomial vol {K + tBi). So 
the magnitude of K should be 

Vol {K) SmijdK) Vi(K) ^ 

8tt 4tt 2 ■ 

In the notation of Theorem [51 the question of validity of the convex magnitude 
conjecture reduces to whether for the solution h G 77® (M®) of the boundary value 
problem (7 — A)®7 = 0 off 77, 7 = 1 on 77, we have 

1 f dAh _ surf (977) Vi(77) 


dAh , „ surf (977) 
— / —;r do = --- 

StT Jqk OV 471 


+ 1 . 


When 77 is the cuboid [—77i,77i] x [—772,772] x [—773 , 773] we have Vi(77) = 2(77i + 
772 + 773 ), so that we expect the magnitude of 77 to be 

R 1 R 2 R 3 {R 1 R 2 + R 2 R 3 + R 3 R 1 ) 


2tt 


+ (77i + R 2 + R 3 ) + 1. 
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In particular in the case of the cuboid do we have 


1 f dAh _ {R 1 R 2 + R 2 R 3 + R 3 R 1 ) 

877 JdK di' 2 tt 


+ (i?i + i?2 + R 3 ) + 1? 


Even when Rj = R for j = 1,2,3, why should the quantity 


— f ^dS 

877 JdK dv 

be a polynomial of degree 2, still less the specific polynomial 3i?^/27r + 3i? + 1? 
One of the difficulties in addressing this problem lies in not having explicit formulae 
for solutions of these types of boundary value problems in the absence of spherical 
symmetry. 


On the other hand, matters may be better if we instead consider ellipsoids. Khavin- 
son and Lundberg note that ellipsoidal harmonics can be used to study elliptic 
boundary value problems in the interior of ellipsoids. It may be possible to use 
them in the current context to obtain explicit solutions of (/ — A)^m = 0 in the 
exterior of an ellipsoid in which might then lead to explicit formulae for the 
magnitudes of ellipsoids in R^. Note however that the mixed volumes, even the 
surface area, of a general ellipsoid in R^ are not expressible as elementary functions 
of the semi-axes, but instead involve elliptic integrals. 


2. The methods presented here can be extended to obtain explicit expressions for 
the magnitude of general compact sets possessing spherical symmetry. 

3. An as-yet-unexplained (empirical) phenomenon is that in the formula for the 

coefficient ao{R) in the n-dimensional setting, with n odd and small, the formula for 
magnitude of the ball of radius R in R”“^ makes a rather mysterious appearance, 
see the previous section. Further observations of this kind have been made by 
Willerton, (private communication), who has also developed a more streamlined 
approach to solving the ODE (fTHll and fitting the boundary conditions. This leads 
to a more symmetric system of linear equations for the unknowns aj. Using the 
new approach to calculating magnitude of Leinster and Meckes Theorem 

4.16), he is able to derive a formula for the magnitude of a ball in odd dimensions 
expressed in terms of the reverse Bessel polynomials. This formula can be re-cast 
in a reasonably succinct and natural form in odd dimensions at least up to and 
including 39 (where the coefficients of the numerator in the magnitude total over 
60, 000 digits). 

4. The methods here can also be used to study the extremal energies or Bessel-like 
capacities 

C™(iL,A) :=inf{||(A/-: h = 1 on K} 

for A > 0 for convex bodies K when to G N and to calculate them explicitly when 
K is a ball. Indeed, a simple scaling argument shows that 

Cm{K,X) = 

and we have shown the asymptotic behaviour of C(^n+i)/ 2 iRi 1) and how to calcu¬ 
late it explicitly for balls above. In particular, C(^n+i)/ 2 {R} as A ^ 0 and 

C'(n+ 1 )/ 2 (A:,A) - X^'^+^'>/‘^Yol{K)/n\LOn as A —>■ 00 . When TO = I and n = 3 for ex¬ 
ample, we can calculate using Theorem[S]and the formula for ipi that Ci{B{0, R), 1) 
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is simply 47r(i?^/3 + i? + 1). So 

C'i(B(0, R),X) = 1) 

= 47rA-^/2(A3/2i?V3 + A^/^i? + 1) 

= 47rAi?V3 + 47ri? + 


5. There remains the possibility that a reformulated version of the convex mag¬ 
nitude conjecture might be true if we replaced magnitude by some closely related 
notion which coincides with it in dimensions one and three, and is strictly smaller 
than it for nontrivial balls in odd dimensions five and above. It would be very 
interesting to find a such a notion. The notion of maximum diversity 

|X|+:=^inf{||/f^(„+„/.(«„) : f € f > 1 on xj 

found for example in HU, HZ] and the references therein, and which is much more in¬ 
tune than magnitude with classical potential theory, unfortunately does not satisfy 
these criteria. One should note the parallel between the definition of maximum 
diversity and certain obstacle problems for higher-order operators (see for example 
[3] and g]) though this does not seem to have been exploited yet.) We thank E. 
Milakis for bringing this parallel to our attention. 


6. In forthcoming work, H. Gimperlein and M. Goffeng [7] use methods from semi- 
classical analysis together with our Theorem [S] to exhibit an asymptotic expansion 
of the magnitude |i?^| ns R ^ oo, where X C K" is open and bounded with 
smooth boundary, when n > 3 is odd. In particular, they prove 


Iflxi = 


nl LJr\ 


, (n + l) Surf(aX)„„_, 

-o—i- ^ 

2 n] uin 


cR 


n-2 


+ 0(i?"-2) , 


where c is proportional to an integral of a curvature of dX. (The coefficient of 
RTi-i differs from the prediction of Leinster and Willerton unless n = 3.) 
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